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ABSTRACT 

We report the first mid-infrared detection of highly disturbed ionized gas in the ultraluminous infrared galaxy 
(ULIRG) IRAS F00183-71 1 L The gas, traced by the 12.81 [Ne ll] and 15.56 [Ne III] lines, spans a ve- 
locity range of -3500 to H-3000 km s"' with respect to systemic velocity. Optical and near-infrared spectroscopic 
studies show no evidence for similarly high velocity gas components in forbidden lines at shorter wavelengths. 
We interpret this as the result of strong extinction (Av=10-50) on the high-velocity gas, which identifies the 
base of the outflow traced in 5007 A [O III] as a plausible origin. Unusual excitation conditions are implied by 
a comparison of the mid-infrared low-excitation neon line emission and the Polycyclic Aromatic Hydrocarbon 
(PAH) emission for a sample of 56 ULIRGs. For IRAS F00183-7111, the neon/PAH ratio is 8 times higher 
than the average ratio. Similar mid-infrared kinematic and excitation characteristics are found for only two 
other ULIRGs in our sample: IRAS 12127-1412NE and IRAS 13451H-1232. Both sources have an elevated 
neon/PAH ratio and exhibit pronounced blue wings in their 15.56 /im [Ne III] line profiles. IRAS 13451H-1232 
even shows a strongly blue shifted and broad 14.32 jjm [Ne v] line. While for IRAS 13451H-1232 the observed 
[Ne lll]/[Ne ll] and [Ne v]/[Ne ll] line ratios indicate exposure of the blue shifted gas to direct radiation from 
the AGN, for IRAS F00183-71 1 1 and IRAS 12127-1412NE the neon fine ratios are consistent with an origin 
in fast shocks in a high-density environment, and with an evolutionary scenario in which strongly blue shifted 
[Ne II] and [Ne III] emission trace the (partial) disruption of the obscuring medium around a buried AGN. The 
detection of strongly blue shifted [Ne v] emission in IRAS 13451H-1232 would then indicate this process to be 
much further advanced in this ULIRG than in IRAS F00183-71 1 1 and IRAS 12127-1412NE, where this Une 
is undetected. 

Subject headings: galaxies: jets and outflows — infrared: ISM — galaxies: ISM — galaxies: active — galax- 
ies: individual (IRAS F00183-71 11, IRAS 12127-1412, IRAS 13451H-1232) 



1. INTRODUCTION 

Large-scale galactic outflows of gas (sometimes called "su- 
perwinds") are seen in galaxies both at low and high red- 
shifts (e.g. Veifleux, Cecil & Bland-Hawthorn 2005; Pettini 
et al. 2002). Superwinds are generated when the kinetic en- 
ergy in the outflows from massive stars and supernovae is 
thermalized, generating a region of very hot (T^^IO^-IO^ K) 
low-density gas in the ISM of a starburst galaxy (Chevalier 
& Clegg 1985). As the bubble breaks out of the disk of the 
galaxy, it will rupture, producing a weakly-collimated bipolar 
outflow into the galaxy halo. As the bubble expands it sweeps 
up and shocks ambient material, creating an outflowing wind 
(Castor, McRay & Weaver 1975; Weaver et al. 1977). Super- 
winds may also be driven, at least in part, by radiation pres- 
sure on dust grains (Murray, Quaetert & Thompson 2005). 
They may also be driven directly by the interaction of jets on 
the dense ISM surrounding an AGN (e.g. Holt, Tadhunter & 
Morganti 2003). 

There is morphological, physical and kinematic evidence 
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for superwinds in nearby starburst and infrared luminous 
galaxies (e.g. Heckman, Armus & Miley 1990). Large-scale 
optical emission-line and associated X-ray nebulae are ubiq- 
uitous in starbursts (Armus, Heckman & Miley 1990; Lehnert, 
Heckman & Weaver 1999; Strickland et al. 2004), and these 
increase in size and luminosity from dwarf starbursts all the 
way to ULIRGs (where they can be tens of Kpc in size — 
(Grimes et al. 2005). In addition, the winds are often seen via 
blue-shifted absorption features (e.g. NaD) against the star- 
burst nuclei (Heckman et al. 2000; Rupke et al. 2005; Martin 
2005) indicating that the neutral, cold ISM is being swept-up 
in the outflow. At high redshifts, superwinds have been seen 
in Lyman break galaxies at z=2-3 (Pettini et al. 2001 ; Shapley 
et al. 2003). 

The ultraluminous infrared galaxy (ULIRG: L(8- 
1000^m)=10'2-i3 Lq) IRAS FOO 183-71 11 is known 
from optical emission line spectroscopy to harbour such a 
galactic-scale outflow (Heckman, Armus & Miley 1990). 
In this paper, we study the mid-infrared evidence for this 
outflow based on the Spitzer-IKS emission line spectroscopy. 
The emission line data is analyzed and presented along with 
new 3.6, 5.8, 70 and 160 /im Spitzer photometric observa- 
tions, obtained to better constrain the infrared spectral energy 
distribution of this interesting ULIRG. 

With an integrated 8-1000 /im luminosity of 8.6±1.3 x lO'^ 
Lq (this work) IRAS FOO 18 3-71 1 1 is one of the most lumi- 
nous ULIRGs discovered by IRAS. Near-infrared imaging by 
Rigopoulou et al. (1999) (Fig. 1) shows a disturbed morphol- 
ogy and what appears to be a single nucleus. Evidence for 
the presence of an accreting central supermassive black hole 
is provided by the radio luminosity, which is in the range of 
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Fig . 1 .— K-band image of IRAS FOOl 83-7 1 1 1 replotted from Rigopoulou 
et al. (1999). North is up. East is to the left. Overlaid on the image are the 
various slits used to obtain spectra of the emission line features discussed in 
the paper. In red and orange is the Spitzer-IRS SH slit in the two nodding 
positions; in blue is the slit used by Heckman, Armus & Miley (1990), in 
green the slit used by Dannerbauer et al. (2005), and in purple the slit used 
by Buchanan et al. (2006) 



powerful radio galaxies and a factor ten in excess of what 
would be expected for starburst galaxies based on the far- 
infrared luminosity (q=1.14; Roy & Norris 1997; Norris et 
al. 1988). The presence of an AGN was recently confirmed 
by Nandra & Iwasawa (2007), who detected a 6.7 keV Fe Ka 
line (Fe XXV) with a large equivalent width, indicative of re- 
flected light from a compton thick AGN. 

The most common star formation tracers, the PAH emis- 
sion bands, were not detected in the 5-12 /im ISO-CAM- 
CVF spectrum (Tran et al. 2001). Instead, the spectrum is 
dominated by a strong silicate absorption band centered at 
9.7 /im, indicative of strong absorption of a background con- 
tinuum source, most likely a deeply buried AGN (Tran et 
al. 2001). One of the PAH emission bands, the C-H bend- 
ing mode at 11.2 /im, was later identified in the Spitzer-IRS 
low-resolution spectrum and was used to constrain the con- 
tribution of mid-infrared detectable star formation to the en- 
ergy budget of IRASF00183-7111 to be <30% (Spoon et 
al. 2004). The same spectrum may also offer indications for 
ongoing buried star formation through the detection of ab- 
sorption bands of crystalline silicate features at 16, 19 and 
23 /im (Spoon et al. 2006). Crystalline silicates are known to 
be forged by young stars in their circumstellar disks and by 
evolved stars in their dense stellar winds (Molster & Kemper 
2005, and references therein). Their presence in the general 
ISM of IRAS F00183-71 1 1 may indicate fierce ongoing star 
formation at mid-infrared optical depths significantly larger 
than one (Spoon et al. 2006). 

Optical long-sHt spectroscopy of IRAS FOOl 8 3-71 1 1 has 
revealed the presence of kinematically highly disturbed ion- 
ized gas extending from 10 kpc West to ~25kpc East of the 
nucleus. This gas, traced in both the 4959 and 5008A [O III] 
lines, is clearly blue shifted between 5 and 20 kpc East and 
shows line widths in excess of 600 km s"' (FWHM) over a 
range of 25 kpc. The double peaked profile of both [O III] 
lines in this range is in agreement with a superwind model in 
which the line emission originates from the thin outer shell 
of shock-heated, swept-up ambient gas. Beyond 25 kpc East, 



TABLE 1 
Photometry of 
IRASF00183-7111 



Band 


Aiest 


Flux density 






(mJy) 


IRAC-1 


2.7 


2.9 ±0.3 


IRAC-3 


4.4 


22 ±2 


IRAS- 12 


9.0 


<60'' 


IRAS-25 


18.8 


133 ± 10" 


IRAS-60 


45 


1200 ± 837" 


MIPS-70 


53 


1500 ± 225 


IRAS- 100 


75 


1190 ± 119" 


MIPS- 160 


121 


540 ± 80 



"IRAS Faint Source Catalog v2 



the gas is much more quiescent, with line widths of ^^200 km 
s"' centered close to systemic velocity (Heckman, Armus & 
Miley 1990). 

Recently, Nandra & Iwasawa (2007) have pointed out that 
the presence of an outflow may also explain the clear eastward 
offset by about 22 kpc of the soft (0.3-l.OkeV) X-ray radia- 
tion detected from IRAS F00183-71 1 1. It also may identify 
the secondary peak at 25 kpc East in the R-band image'^ by 
Drake et al. (2004) to be related to this outflow (Nandra & 
Iwasawa 2007). A corresponding K-band peak is absent in 
the 2.2 /im image (rest frame 1.65 /im) by Rigopoulou et al. 
(1999) (Fig.l). 

In this paper we report the first identification in a ULIRG of 
broad, asymmetric and somewhat blue shifted emission lines 
of [Ne II] at 12.81 /im and [Ne III] at 15.56 /tm in the high- 
resolution Spitzer-\RS specti'um of IRAS F00183-71 1 1. We 
compare our detection to available optical and near-infrared 
spectroscopy to investigate the processes responsible for this 
broadening. Our paper is organized as follows. In Section 2 
we present our observations of IRAS F00183-71 1 1. In Sec- 
tion 3 we present the 1-120 /tm SED and an analysis of the 
mid-infrared emission line spectrum. In Section 4 we analyze 
line profiles of two other ULIRGs which display strong blue 
shifts in their fine-structure neon line profiles. Results are dis- 
cussed in Section 5 and conclusions are provided in Section 6. 
Additional results for IRAS F00183-7111 are presented in the 
Appendix. Throughout this paper we assume Hq = 7 1 km s"' 
Mpc-i, = 0.27, Oa = 0.73 and ilg = 0. 

2. OBSERVATIONS AND DATA REDUCTION 

2. 1 . IRAC and MIPS photometry 

Spitzer Space Telescope (Werner et al. 2004) IRAC (Fazio 
et al. 2004) 3.6 and 5.8 pm photometry was obtained on 2006 
August 14 (AOR key 17517824). IRAS FOOl 8 3-71 1 1 was 
imaged at 5 dither positions using 2 sec ramp times. The 
source flux density was measured by integrating over a 12" 
radius aperture placed upon the SSC pipeline-produced (S14) 
post-BCD mosaic. We estimate the uncertainty in the pho- 
tometry due to the background to be ^2%. The absolute cali- 
bration uncertainties of the IRAC post-BCD images are 10%. 
As part of the same program we obtained MIPS (Rieke et al. 
2004) 70 and 160 /im photometry on 2006 June 17 (AOR key 
17528320). IRASF00183-7111 was imaged for 1 cycle at 
70 /im and for 4 cycles at 160 /im using 3 sec ramps. The 

^ at the redshift of the source the R-band filter includes the 5008A [O III] 
line (Drake et al. 2004) 
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Fig. 2. — The near-to-far-infrared spectral energy distribution of 
IRAS FOO 183-7111. The 4-29 ^tm SpiU,er-VAS low-resolution spectrum 
(Spoon et al. 2004) is shown as a continuous line. Photometric points are 
labeled with the following symbols: circles: ESO-NTT J & K (Rigopoulou 
et al. 1999), squares: IRAC 3.6 & 5.8 /im, triangles: IRAS 60 & 100 ^tm, 
diamonds: MIPS 70 & 160 /im. Error bars on the photometric points are 
too small to be visible. Inset: detailed view of the spectral range around 
the 9.7 \xm silicate absorption profile. Note the emission lines of H2 S(3) at 
9.66 \xm and [Ne 11] at 12.81 ixm. 

source flux densities were measured from the SSC pipeline- 
produced (SI 4) unfiltered post-BCD images by integrating 
the flux within 8 and 20" radius apertures at 70 and 160 \im, 
respectively. Aperture corrections were applied to recover the 
fuU point-source flux densities. We estimate the uncertainty 
on the 70 and 160/xm flux densities due to the background 
to be around 1 and 2%, respectively. The absolute calibra- 
tion uncertainties of the MIPS post-BCD images are 10-20%. 
We present the results of our photometric observations in Ta- 
ble 1. The calculation of infrared luminosities which utilize 
this photometry are described in Section 3.1. 

2.2. IRS spectroscopy 

High resolution (Av ~500 km s"'; R~600) 10-37 /zm IRS 
(Houck et al. 2004) spectra of IRAS F00183-711 1 were ob- 
tained at two epochs: 2005 June 2 (AOR key 6651904) and 
2008 June 1 (AOR keys 25289472 and 25289984). In botii 
cases the source was acquired using a high-accuracy blue 
peak-up and observed in staring mode. The 10-19.5 /zm por- 
tion of the spectrum was obtained in 9 cycles of 120 sec 
integrations in IRS Short High (SH; slit size 4.7" x 11") at 
R\=91°. The 19. 3-37.0 /xm portion of the spectrum was ob- 
tained in 2 sets of 35 cycles of 60 sec integrations, alternated 
with sky background observations (30 cycles of 60 sec), in 
IRS Long High (LH; slit size ll.l"x22.3") at PA=-15°. Fig- 
ure 1 shows the SH slit overplotted on the K-band image of 
the source (Rigopoulou et al. 1999). The LH slit is not shown 
as it is more than twice as wide and long as the SH sUt. 

Data reduction started from 2D droop level (i.e. pre- 
flatfield) images provided by the S 17.2 pipeline of the Spitzer 
Science Center (SSC). For SH, in a first step we used 
IRS CLEAN to interpolate over the bad and 'rogue' pixels. For 
LH, this step was preceeded by subtraction of the contempo- 
raneous sky background images. We then used the SMART 
data reduction package (Higdon et al. 2004) to extract the ID 
target spectrum in 'full slit' mode. The same procedure was 
followed for a large number of observations of the calibration 
star ^ Dra. The final spectrum was obtained by multiplying 
the spectrum of IRASF00183-7111 by the relative spectral 
response function created from the ratio of the observed spec- 
trum of ^ Dra to the stellar reference spectrum (Cohen et al. 
2003, G.C. Sloan priv.connm.). Note that given the absence of 



TABLE 2 

Synthetic photometry for 
IRAS FOG I 83-71 I 1 



Band 


Arest 


Flux density 


Uncertainty 




(Aim) 


(mjy) 




IRAC-1 


3.6 


8.6 


1% 


IRAC-2 


4.5 


29.6 


1% 


IRAC-3 


5.8 


47.3 


1% 


IRAC-4 


8.0 


48.6 


1% 


MIPS-24 


24 


370 


1% 


IRAS-25 


25 


360 


1% 


IRAS-60 


60 


1410 


1% 


MIPS-70 


70 


1190 


2% 


IRAS- 100 


100 


890 


3% 


MIPS- 160 


160 


490 


5% 


SCUBA-450 


450 


20 


6% 


SCUBA-850 


850 


2.0 


6% 



Notes. — Rest-frame synthetic fluxes measured 
from the model spectrum of IRAS FOO 1 83-7 1 1 1 (see 
Section 3.1). The statistical uncertainty in the syn- 
thetic photometry is computed from the uncertainty 
in the total fitted flux 



contemporaneous SH background observations, the SH por- 
tion of the spectrum is not background-subtracted. 

The wavelength calibration of our spectrum is accurate to 
about 1/5 (SSC, priv. conrai.) of a resolution element (500 
km s~^), which amounts to 100 km s"^ . 

3. ANALYSIS 
3.1. Infrared SED 

With the addition of the IRAC 3.6 /im and MIPS 70 and 
160 /xm photometry it is possible to fit the infrared SED of 
IRAS FOO 183-71 11 over the entire rest frame wavelength 
range of 1-120 /zm. Especially important is the addition of the 
MIPS 160 /im point, which extends the sampling of the rest 
frame SED from 75 /im (pie-Spitzer) to 120 /im (see Fig. 2), 
allowing a more rehable extrapolation of the SED into the 
submilUmetre range. For the latter, we have used the multi- 
component spectral energy distribution decomposition code 
ofMarshaUetal. (2007). 

Based on our fit, we have computed the following integrated 
rest frame infrared luminosities for IRAS F00183-71 1 1: L(l- 
1000/xm)=9.9±1.3xl0i2LQ, L(8-1000/im)=8.6±1.3xl0i2 
Lq, L(40-500/im)=4.4±0.8xl0i2 l(4o_ 

120/im)=4.0±0.6xl0i^ Lq. We have further used the 
results of the fit to derive rest frame filter flux densities for 
the IRAC, MIPS, IRAS and SCUBA bands. These values 
are tabulated in Table 2. Given the 15% absolute calibration 
uncertainty of the MIPS 70 and 160 /xm data, and the fact 
that the MIPS 160 /xm data is the longest wavelength point 
constraining our fit, there is clearly some uncertainty in our 
extrapolation of the flux density out to 1000 ^m. The spectral 
shape of the extrapolated flux density beyond 160/im is a 
function of the chosen power-law spectral index of the far-IR 
dust opacity. The above stated uncertainties in the integrated 
luminosities incorporate both the systematic uncertainty 
resulting from our choice of this dust opacity spectral index 
as well as the statistical uncertainty resulting from the MIPS 
calibration uncertainty. 

The ratio of the observed IRAS-25 to IRAS-60 filter fluxes 
is commonly used to classify ULIRGs as "warm" (25/60 /xm 
> 0.2) or "cold" (25/60 /xm < 0.2), separating ULIRGs be- 
lieved to be AGN-dominated from those which are thought 
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Fig. 3. — Comparison of the emission line profiles for the six detected emission lines in the spectrum of IRAS F00183-71 11. The left column shows the S(l), 
S(2) and S(3) pure rotational H2 lines, the right panel the [Ne II], [Ne III] and [S III] lines. All six lines have been observed at the same spectral resolution S=600. 
The [Ne II] and [Ne III] line profiles are clearly asymmetric and broadened beyond the instrumental spectral profile (FWHM = 500 km s"'). Their profiles are fit 
by two Gaussian components (red and blue dotted lines, respectively), the sum of which is shown as a red or blue continuous line 



to be powered predominantly by star formation (de Grijp et 
al. 1985; Sanders et al. 1988). For the classification to be 
meaningful also for higher redshift targets like IRAS F00183- 
7111, rest frame synthetic IRAS fluxes should be used. For 
IRASF00183-7111, this results in an IRAS-25 to IRAS-60 
flux ratio of 0.26, which places it among the class of "warm", 
AGN-dominated ULIRGs, in agreement with results obtained 
at other wavelengths (see Section 1). 

3.2. Emission lines 

The high-resolution Spitzer-W.Si spectrum of 
IRAS FOO 18 3-71 1 1 covers the rest frame wavelength 
range of 7.5-27.8 jim. In Table 3 we list the line fluxes and 
3(7 upper limits for 17 emission lines, as well as for the 
1 1.2 jim. PAH emission feature. For the latter, the feature flux 
was measured above a spline continuum defined shortward 
of 11.1 /im and longward of 11.6 /im. The only lines that 
are detected in our spectrum are the S(l), S(2) and S(3) pure 
rotational lines of molecular hydrogen, and the 12.81 /im 
[Ne 11], 15.56 /im [Ne ill] and 18.71 /im [S III] fine-structure 
lines. In Fig. 3 we present their continuum-subtracted line 



profiles as a function of velocity relative to systemic^. 

Of the six detected emission lines, the two neon fine- 
structure lines appear clearly resolved beyond the instrumen- 
tal resolution (500 km s"' ; R=6QQi), most notably the [Nell] 
line, at an observed FWHM of 1500±100km s"' (Table 3). 
The other emission lines have observed FWHM ranging from 
490±50 to 950±130km s"'. In what follows any FWHM 
values quoted in the text refer to the true FWHM, i.e. the 
observed FWHM corrected for the instrumental profile. 

3.3. Neon line profile analysis 

Extracting a [Ne 11] line profile requires careful removal of 
the underlying spectral structure. As illustrated in Fig. 4 by 
the starburst galaxy NGC2146 and the ULIRG NGC6240, 
the [Nell] line is superimposed on the 12.7 /im PAH emis- 
sion feature, which arises over a velocity range of -13500 to 

' In the absence of CO or 21cm redshift determinations, we have used the 
mid-infrared S(l) and S(3) molecular hydrogen lines to determine the reces- 
sion velocity of IRAS FOO 18 3-7 111 to be 98400±200km s"' (z=0.328 ± 
0.001). This redshift is 360km s"' higher than the value adopted by Heck- 
man, Armus & Miley (1990) and 180 km s"' higher than the value used by 
Roy & Norris (1997) 
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TABLE 3 

Emission FEATURES in the spectrum of IRAS FOO 183-71 1 1 



Feature ID 


^rest 


; 

Flux 


ODS. rWriM 


tl Lie r W rIM 






(10^^' Wcm"^) 


(km s ^ ) 


(km s ^) 


[Ne VI] 


7.65 


<0.25 






Ht S(4) 


8.03 


<0.26 






[Ar III] 


8 99 


<0 07 






Ht S(3) 


9.66 


n ^9 -1- 0^ 


1 1. )\ 1 


410di50 


rs TVl 

[O 1 V J 


10.51 


<0.08 






PAH 


1 1.2 


2.7 it 0.4 






Ht S(2) 


12.28 


0.59 ± 0.1 


950±130 


810ibl20 


HuQ 


12.37 


<0.12 






[Ne II] 


12.81 


6.3 ± 0.3 


1500±100 


1410±100 


[Ne V] 


14.32 


<0.13 






[CI II] 


14.37 


<0.14 






[Ne III] 


15.56 


2.1 ± 0.2 


1300±120 


1200±120 


Hi S(l) 


17.04 


0.83 ± 0.05 


660±50 


430±55 


[Fe IT] 


17.94 


<0.10 






[S III] 


18.71 


0.23 ± 0.03 




<ZUU 


[Ne V] 


24.30 


<0.20 






[OIV] 


25.89 


<0.35 






[Fe n] 


25.99 


<0.31 






[OlII]" 


0.5007 


0.45 




890 


[0 III]'^ 


0.5007 






600-1000 


Ha'' 


0.6565 


0.51 ±0.03 




720 


[Fe 11]" 


1.257 


0.228: 




900: 




1.282 


0.115: 




600: 


[Fe 11]" 


1.644 


0.376: 




1000: 


"Dannerbauer et al. (2005) 






''Buchanan et al. (2006) 






'^Heckman, Armus & 


Miley (1990) 






Notes. — 


- Line fluxes and 3cr upper limits above the separator line 



are measured from our Spitzer-IRS spectrum, those below have been taken 
from the literature. The uncertainties in our line fluxes are the larger of 
either the error in the Gaussian fit or half the difference of the two nodding 
positions. The FWHM in column four (five) is from a Gaussian fit to the 
data, before (after) removal of the instrumental profile, in quadrature. 



+3500 km s ' relative to the [Ne ll] line. However, unlike 
NGC2146 and NGC6240, the spectrum of IRASF00183- 
7111 shows hardly any PAH emission, with the exception of 
the PAH feature at 11.2/im, which is detected at very low 
equivalent width. We conclude that the remarkable width of 
the 12.81 /im [Ne ll] line can not be ascribed to inproper sub- 
traction of the underlying spectral structure. In fact, the re- 
sulting profile is consistent with that of the 15.56 /im [Ne III] 
line. 

The shape of the [Ne ll] line profile is asymmetric, with 
the peak emission arising very close to systemic velocity. To 
characterize the asymmetry we have decomposed the line pro- 
file into two Gaussian components (Fig. 3). The resulting fit 
is good. We find the dominant component to have a FWHM 
of 2570±90kms-i, centered at a blueshift of 400±50km s"' 
w.r.t. systemic. The narrow component has a FWHM of 
610±60kms"' and a redshift of 60±20kms"' w.r.t. sys- 
temic. The fit components are shown overplotted in red onto 
the observed [Ne ll] profile in the top-right panel of Fig. 3. 
Further fit details are provided in Table 4. 

Unlike the 12.81 fim [Ne ll] Une, the 15.56 [Ne III] Une 
resides in a relatively uncomplicated part of the spectrum. 
We determined its profile by defining and subtracting a local 
third-order polynomial continuum. We then performed a 2- 
component decomposition in the same way as for the [Ne ll] 
line. The narrow component has a FWHM of 460±90km s"' 
and is centered at the systemic velocity, whereas the broad 
component has a FWHM of 2850±240km s"' and is blue 
shifted by 450±100km s"'. The fraction of the [Ne III] flux 
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Fig. 4. — Comparison of the continuum-subtracted spectral structure 
around the 12.81 ^m [Ne II] line for four galactic nuclei, at a spectral 
resolution S=600. The spectral features included are the S(2) line at 
12.28 /im, the 12.7/im PAH emission feature (extending from -13,500 to 
+3,500 km/s) and the [Ne II] line. The latter is clearly broadened in the 
spectra of IRAS FOO 183-71 1 1 and NGC6240. Note that the 12.7 /tm PAH 
emission feature appears to be absent in the spectra of IRAS F00183-71 1 1 
and IRAS 08572+3915NW 

emitted in the broad component is 3/4, slightly lower than the 
fraction found for the [Ne ll] line, 4/5. As a result, the [Ne 
lll]/[Ne II] ratios are 0.31 ±0.03 for the broad component and 
0.39±0.08 for the narrow component. 

3.3.1. Comparison to other line measurements for 
IRAS FOOl 83-71 11 

The most detailed spectral line study that we are aware of 
focusses on the [O III] doublet at 4959 and 5007 A and was 
published by Heckman, Armus & Miley (1990). They ob- 
served IRASF00183-7111 in a 2"xl00" slit along PA 90° 
(see Fig. 1) at a spectral resolution R=\9QQ. In Fig. 5 we 
compare the 5007 A [O III] profile obtained at a position 3" 
(13kpc) East of the nucleus to that of our galaxy-integrated 
spectrum, after degrading the spectral resolution of the [O III] 
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Fig. 5. — Comparison of the galaxy-integrated line profile of 12.81 fim 
[Ne II] (red) and 5007 A [O III] (green) at a position 3" East of the nucleus 
of IRAS FOO 183-71 11, where the 0^+ gas reaches its highest blue shift. The 
spectral resolution of the double-peaked [O III] profile (R=1900) has been de- 
graded to the resolution of the 5prYzer-IRS-SH measurement (R=600). Over- 
plotted are the 2-component Gauss fits to the [Ne II] and [Ne III] line profiles. 
The individual components are shown as dashed lines, their sums as dotted 
lines 

line from 1900 to 600. According to Fig. 12 of Heckman, Ar- 
mus & Miley (1990), at this position the centroid of the [O 
III] line profile has its highest blueshift (700km s"' ; using our 
preferred redshift of 0.328±0.001) at a FWHM -600 km s"', 
typical for the entire region extending from lOkpc West to 
20kpc East of the nucleus. Clearly, the [Ne ll] and [Ne III] 
emission extends to much higher velocities than the [O III] 
line emission does. 

More recently, Buchanan et al. (2006) published line fluxes 
and FWHMs for 11 optical lines observed at /?=1000 in a 
2"x9" slit oriented along PA=329° (C.L. Buchanan, private 
communication). For the Ha line, which is strongly blended 
with the [ N ll] line on either side, they infer a FWHM of 720 
km s~' (compared to 700 km s"' reported by Armus, Heckman 
& Miley 1989). For the 4959/5007 A [O III] doublet, they re- 
port FWHMs of 800-900 km s"' (Table 1), which are com- 
parable to those given by Heckman, Armus & Miley (1990) 
in their Figure 12, even though the slit orientations in the two 
studies in question differ by 239 degrees (Fig. 1). 

We also compared our mid-infrared results to near-infrared 
H and K band spectroscopy by Dannerbauer et al. (2005) ob- 
tained using a l"x290" sht oriented at PA=90° (Fig. 1). Their 
data have similar spectral resolution as ours. The two detected 
[Fe II] lines at 1 .257 and 1 .644 /im have FWHMs in the range 
900-1000 km s"', while the 1.282 ^m Pa/3 Hne appears to be 
unresolved. Again, like in the optical range, none of the lines 
have FWHMs as large as the mid-infrared neon lines (see Ta- 
ble 3). 

3.4. AGN tracing mid-infrared lines 

The rest wavelength range covered by the Spitzer IRS high- 
resolution modules includes three high-ionization lines com- 
monly used as tracers of AGN activity: 7.65 fim [Ne VI] and 
14.32 &24.32/xm [Ne v] (e.g. Sturm et al. 2002). None of 



them are detected (Table 3). The 3cr limit on the fraction of 
8-1000 /im infrared power contained in the 14.32 /im [Ne v] 
line is L[Nev]/LiR < 1-3 x 10"^, which is 1-2 orders of magni- 
tude smaller than the detections for classical AGNs, but sim- 
ilar to the limits derived for other ULIRGs with deep 9.7 /im 
silicate features like IRAS 08572H-3915, IRAS 15250H-3609 
and IRAS 20100^156, and for the type-1 ULIRG/LoBAL- 
quasar Mrk231 (Armus et al. 2007). For the limit on 7.65 /im 
[Ne VI] no ULIRG comparison data exists. However, the 3a 
upper limit L[Nevi]/LiR < 2.6 x 10"^ is at least an order of mag- 
nitude below what is found for classical AGNs (Sturm et al. 
2002). 

3.5. PAH emission 

The quality of our Spitzer IRS high-resolution spectrum al- 
lows us to measure the flux in the 11.2 fim PAH feature more 
accurately than has been possible from the low-resolution 
spectrum (Spoon et al. 2004). The latter showed strong dis- 
agreements between the spectra obtained in nod positions 1 
and 2 at the short wavelength end of IRS Long-Low order 2. 
The new feature flux, measured from the IRS SH spectrum 
is tabulated in Table 3 and is 25% smaller than previously 
reported by Spoon et al. (2004) based on the low-resolution 
spectrum. We also remeasured the upper limit on the 6.2 fim 
PAH feature from a re-extracted low-resolution spectrum, us- 
ing the same continuum pivots at 5.65, 6.07 and 6.50 fim to 
define the local spline continuum as in Spoon et al. (2004). 
Given the clear reduction in noise and artefacts between IRS 
pipeline versions 9.1 and 14.2, the new 3a upper limit on the 
6.2 /im PAH flux is 1.9x 10"^' W cm"-, which is a factor 
of 2 lower than previously reported by Spoon et al. (2004). 
Note that we do not attempt to measure the flux in the 7.7 fim 
PAH band, as this emission feature — if present — would be 
strongly blended with the silicate feature. 

Based on the new PAH fluxes we compute an upper limit 
on the ratio of 6.2-to-l 1.2 fim PAH of 0.70. This is at the low 
end of the range found for starburst nuclei and ULIRGs in our 
sample {2a below the average ratio of 1 .4), but not untypical 
for a sample of IR-faint LINERs (Lir=10^-^^-'°-5^) studied by 
Sturm et al. (2006), and for the spiral arms of M 5 1 and the 
halo of M 82 (GalHano et al. 2008). 

The low 6.2-to-l 1.2 /im PAH ratio may indicate a high ra- 
tio of neutral to ionized PAHs, which would imply a low 
ratio of the ionizing radiation field to the electron den- 
sity, and physical conditions similar to those in the Galac- 
tic reflection nebula NGC2023 (Gafliano et al. 2008). In 
IRAS FOO 18 3-71 1 1 these conditions could plausibly exist 
outside the deeply buried nucleus. Alternatively, the low ob- 
served 6.2-to- 1 1 .2 /im PAH ratio for IRAS FOOl 83-7 1 1 1 may 
indicate destruction of small PAHs exposed to hard radiation 
from an AGN. This would selectively lower the fluxes emit- 
ted in the PAH emission bands at 3.3, 6.2 and 7.7 fim relative 
to those at 11.2/im and beyond (Smith et al. 2007). Given, 
however, the absence of detections of high ionization lines, 
indicative of a hard AGN radiation field, this scenario may 
not hold for IRAS F00183-71 1 1. 

Note that strong extinction on the PAH emitting region can 
be ruled out as a third explanation for the low observed 6.2-to- 
11.2 fim PAH ratio, as this would have the opposite effect on 
the PAH ratio: the attenuation on the 1 1.2 fim PAH feature is 
stronger than the attenuation on the 6.2 fim PAH feature (e.g. 
Chiar & Tielens 2006). 

4. OTHER ULIRGS SHOWING RESOLVED MID-IR NEON LINES 
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Fig. 6. — Comparison of the galaxy-integrated line profiles of 12.81 fim 
[Ne II] (red) and 15.56 ^im [Ne III] (blue) for IRAS 12127-1412NE at a spec- 
tral resolution of S=600. Overplotted are 2-component Gauss fits to the line 
profiles. The individual components are shown as dashed lines, their sums as 
dotted lines 

The discovery of high velocity [Ne ll] emission in 
IRAS FOO 183-71 1 1 raises the question whether there are 
other luminous infrared galaxies that show evidence for dis- 
turbed (FWHM>800 km s"') neon gas in their Spitzer-IRS 
spectra. If so, do these galaxies have properties similar to 
IRAS FOO 18 3-711 1? To address this question we have in- 
spected the line profiles of the low-ionization emission lines 
12.81 Aim [Ne ll] and 15.56 ^im [Ne III] of the 53 ULIRGs in 
the sample of Farrah et al. (2007), as well as the archival spec- 
tra of the ULIRG IRAS 12127-1412NE and the HyLIRGs 
IRAS F09104H-4 109 and IRAS F15307H-3252. In our seai'ch 
we distinguish two classes. One consisting of sources with 
symmetric resolved neon line profiles, and another with 
sources which show asymmetric resolved line profiles. 

The class of asymmetric resolved [Ne ll] and [Ne III] line 
profiles consists of only two sources besides IRASF00183- 
7111. In all three cases the line asymmetry takes the form of 
a blue wing. 

4.1. IRAS 12127-1412NE 

One source is IRAS 12127-1412NE (Fig. 6), the brighter of 
the two nuclei in this widely separated LINER system (pro- 
jected nuclear separation 9", or 21.9 kpc; Veilleux, Kim & 
Sanders 1999). In the near-infrared, Imanishi, Dudley & Mal- 
oney (2006) classified the NE nucleus as a buried AGN based 
on the absence of 3.3 /im PAH emission and the presence of 
3.0 iim water ice and 3.4 jjm hydrocarbon absorption bands. 
The deeply buried nature of the nuclear power source(s) is 
confirmed by the presence of a deep 10 /im silicate feature 
in the Spitzer-\RS> low-resolution spectrum (Imanishi et al. 
2007). 

On 2008 July 10, we reobserved the source in staring mode 
in IRS-SH, using an integration time of 12 cycles of 120 sec. 
Sky-subtraction was performed using a contemporaneous sky 
background observation of 6 cycles of 120 sec. Data reduction 
proceeded along the same steps as detailed in Sect. 2.1. 

Like for IRAS FOO 18 3-71 1 1, the asymmetry of the neon 




Fig. 7. — Comparison of the galaxy-integrated line profiles of 12.81 ^m 
[Ne 11] (red), 15.56 ^tm [Ne III] (hlue) and 14.32 ^tm [Ne V] (black) for 
IRAS 13451+1232 at a spectral resolution of R=600. Overplotted are 2- 
component Gauss fits to the neon line profiles. The individual components 
are shown as dashed lines, their sums as dotted lines. The [Ne V] line has 
been scaled to match the blue wing of the [Ne III] line profile. Also shown 
is the sum of the 3-component Gauss fit to the 5007 A [O III] Une (green) as 
presented by Holt, Tadhunter & Morganti (2003), degraded to a resolution of 
i?=600, and the Paa profile (purple) obtained by Veilleux, Sanders & Kim 
(1997)atR=650 

emission lines of IRAS 12127-1412NE may be character- 
ized using a 2-component gaussian fit. For the [Ne ll] line 
the narrow component is practically unresolved (FWHM = 
175±45km s"') and is centered close to systemic velocity"' 
(-80±10km s"'). The other component has a FWHM of 
1670± 100km s"' and is centered at a blue shift of 530±60km 
s"' (Fig. 6). The fraction of the [Ne ll] line flux emitted in 
the broad component is slightly smaller than that in the nar- 
row component. This is less than for IRAS F00183-71 1 1, for 
which about 4/5 of the line flux is associated with the broad 
component. We also examined the 15.56 /im [Ne III] pro- 
file. Its shape is not as well constrained as that of the [Ne ll] 
line due to a five times lower line flux, but the profile clearly 
shows the presence of a blue wing. A 2-component Gauss fit 
results in an unresolved narrow component centered at sys- 
temic velocity (-7±40km s"') and a broad component with 
FWHM of 1500±340km s"' blueshifted by 650±200km s"' . 
We find that the contribution of the broad component to the 
[Ne III] line profile is far stronger than to the [Ne ll] pro- 
file, 3/4 as compared to 1/2 (Fig. 6). This implies a dif- 
ference in excitation of the narrow and broad components. 
For the narrow component we find a [Ne lll]/[Ne ll] ratio of 
0.095±0.04, while for the broad component we compute a ra- 
tio of 0.29±0.08 Unfortunately, a comparison to the line pro- 
file of the optical [O III] doublet is not possible as the equiva- 
lent widths of these lines in the spectrum of Veilleux, Kim & 
Sanders (1999) are too low for a meaningful analysis of their 
profiles. 

We infer the systemic velocity of IRAS 12127-1412NE from the 
line centers of the H2 S(l) and S(3) lines and find 39870±200km s"' 
(z=0.1330±0.0004), in good agreement with the redshift of 0.133 quoted 
by Veilleux, Kim & Sanders (2002) 
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TABLE 4 

Fit results for mid-infrared neon emission lines 



Feature 


Component 


Ac 


Vc 


obs. FWHM 


true FWHM 


Flux 






(/im) 


(kms-i) 


(km s-i) 


(kms"') 


(10-2' Wcm-2) 




IRAS 00183-7111 


LNe IlJ 


n 


12.816 


+ 60 ± 20 


790 ± 60 


610 ± 60 


1.4 ±0.16 


[Ne II] 


b 


12.796 


^400 ± 50 


2620 ± 90 


2570 ± 90 


4.9 ± 0.20 


[Ne III] 


n 


15.555 


± 30 


680 ± 90 


460 ± 80 


0.55± 0.09 


[Ne III] 


b 


15.532 


-450 ± 100 


2890 ± 240 


2850 ± 240 


1.5 ±0.12 


[Ne V] 




14.322 








<0.13 








IRAS 12127 


-1412Nli 






LNe llj 


n 


12.818 


-80 ± 10 


530 ± 20 


175 ± 45 


1.68 ±0.08 


[Ne II] 


b 


12.799 


-530 ± 60 


1740 ± 100 


1670 ± 100 


1.47 ±0.10 


[Ne III] 


n 


15.565 


-7 ± 40 


430 ± 120 


<200 


0.16 ±0.07 


[Ne III] 


b 


15.531 


-650 ± 200 


1580 ± 360 


1500 ± 340 


0.42 ±0.10 


[Ne V] 




14.322 








<0.08 


IRAS 13451+1232 


[Ne IlJ 


n 


12.808 


-120 ± 10 


790 ± 40 


610 ± 50 


4.3 ± 0.30 


[Ne II] 


b 


12.758 


-1290 ± 215 


2520 ± 350 


2470 ± 340 


2.5 ± 0.40 


[Ne III] 


n 


15.545 


-200 ± 25 


840 ± 70 


670 ± 70 


1.9 ± 0.25 


[Ne III] 


b 


15.506 


-950 ± 85 


2410 ± 110 


2360 ± 110 


4.1 ±0.30 


[Ne V] 


h 


14.263 


- 1 240 ± 80 


2360 ± 190 


2300 ± 190 


1 .9 ± 0. 1 3 



Notes . — Results from a 2-component gaussian decomposition of the 12.81, 14.32 and 15.56 /tm 
neon emission line profiles into narrow (n) and broad (b) components. Ac denotes the central wave- 
length of the component and corresponds to a velocity shift Vc of the component with respect to 
systemic velocity. The formal error on the velocity shift does not include corrections for the uncer- 
tainty in the wavelength calibration, which is 1/5 of a resolution element (100 km s"'). The FWHM 
in column five (six) is measured before (after) removal of the instrumental profile (500 ± 60 km s"' ; 
Dasyra et al. 2008), in quadrature 



4.2. IRAS 13451+1232 

The third ULIRG with asymmetric neon line profiles is 
IRAS 13451+1232", also known as PKS 1345+12 (Fig. 7). 
The source, which is optically classified as Seyfert-2 
(Veilleux, Sanders & Kim 1997), comprises two nuclei sep- 
arated by 1.8" (4.3 kpc; unresolved at the spatial resolution 
of Spitzer-IRS). Both the low resolution near-infrared 2.6- 
3.65 /im (Imanishi, Dudley & Maloney 2006) and the Spitzer 
mid-infrared spectrum is virtually featureless, characteristic 
of AGN heated hot dust. At the higher resolving power of 
IRS-SH, emission Une profiles of 12.81 nm [Ne ii], 15.56 /xm 
[Ne III] and 14.32 /xm [Ne v] can be discerned on the contin- 
uum, the clearest of which is the is the [Ne III] line. We fit 
its profile with a broad component of FWH]VI 2360±1 10km 
s"' blue shifted by 950±85km s"' and a narrow compo- 
nent of FWHM 670±70km s"' blueshifted by 200±25km 
s^^ The broad component is responsible for 2/3 of the 
[Ne III] emission, significantly more than for IRAS F00183- 
7111 but somewhat less than for IRAS 12127- 1412NE. For 
the [Ne ii] Une profile the 2-component decomposition re- 
sults in a narrow component with a FWHM of 610±50km 
s"^ blueshifted by 120±10km s"^ and a broad component 
with FWHM 2470±340km s"' blueshifted by 1290±220km 
s"'. The broad component carries 1/3 of the total [Ne ll] 
line flux, significantly less than for IRAS FOO 18 3-71 1 1 and 
IRAS 12127-1412NE. The third line, the [Ne V] coronal line 
at 14.32 /Lim, only has a blue shifted component, centered at 
a -1240±80km s~\ which we fit with a Gaussian profile of 

" We adopt z=0.12174±0.00001 as the systemic velocity for 
IRAS 13451+1232. This velocity corresponds to the velocity center of the 
narrow component in the 3-component fit to the 4959 & 5007 A [O III] dou- 
blet in the R=12()0 spectrum of Holt, Tadhunter & Morganti (2003). Using 
this redshift, the line centers of the H2 S(l) and S(3) lines occur at 40±30 
and 20±35 km s ' , respectively. Note that our adopted redshift is 200 km s"' 
lower than the value derived from CO observations (Evans et al. 1999) 



FWHM 2300±190km s*'. The presence of [Ne V] emission 
only at blueshifted wavelengths is consistent with the differ- 
ence in excitation between the narrow and broad components 
as measured from the [Ne iii]/[Ne ii] ratio. For the narrow 
component this ratio is 0.44±0.07, while for the broad com- 
ponent it is 1.6±0.3. The [Ne V]/[Ne II] ratio for the broad 
component is 0.76±0. 14. 

In Fig. 7 we also show the sum of a 3-component fit to 
the profile of the 5007 A [O III| line, centered on the nu- 
cleus (Holt, Tadhunter & Morganti 2003). The spectral res- 
olution was degraded from /?=1200 to 600 for a proper com- 
parison to our IRS-SH data. As a result of the degrada- 
tion, the narrowest (FWHM=340±23 km s"') and weakest 
of the three kinematic components, centered at systemic ve- 
locity, gets blended with the intermediate width component 
(FWHM=1255±12km s"') centered at -402±9km s"'. The 
broadest, most blueshifted component (FWHM=1944±65km 
s"'), centered at -1980±36km s~\ can however still be rec- 
ognized as a separate component. We also overplot the line 
profile of Pact (Veilleux, Sanders & Kim 1997) with a narrow 
component of FWHM=773 km s"' and a broad component of 
FWHM=2588kms-i. The overlay suggests that the 15.56 /itm 
[Ne III] line traces more or less the same blue shifted kine- 
matic components as the 5007 A [O iii] and Paa lines. The 
comparison further shows that, unlike IRAS F00183-71 1 1, 
there is no evidence for the 5007 A [O ill] line to cover a 
more limited velocity range than the 12.81 fim [Ne ll] and 
15.56 fim [Ne iii] lines (Figs. 5 & 7). The latter is in agree- 
ment with the differences in the distribution of the obscuring 
medium as probed by the 9.7 fxm silicate absorption features 
(see Section 5.4). 

5. DISCUSSION 
5.1. AGN activity in IRAS FOO 183-71 11 



9 



Q- 1.0 
E 



13451. 



^FOOIE 



i(rk463 





Silicote Strength 
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of 2.2. The horizontal dotted line indicates the average neon/PAH ratio for low-excitation ULIRGs. IRAS F00183-711 1 is more than a factor 8 above this value, 
indicating that it is overluminous in its neon lines compared to its PAH emission. The same is true to a lesser extent for the ULIRGs IRAS 12127-1412NE and 
IRAS 13451+1232. Galaxy types are labeled as follows: ULIRGs: red, starburst galaxies: blue, Seyfert-2 galaxies: dark green, Seyfert-1 galaxies and QSOs: 
light green 



Observations at various wavelengths suggest that 
IRAS FOO 18 3-71 11 derives a substantial fraction of its 
bolometric power from accretion onto a central black hole. 
Especially important in this respect is the detection by 
Nandra & Iwasawa (2007) of a 6.7 keV Fe Ka line (Fe XXV) 
with a large equivalent width, indicating the presence of a 
Compton-thick AGN in IRAS F00183-71 1 1. Another clear 
indication is a higher than expected radio power based on 
the radio-far-infrared correlation for starburst galaxies (Roy 
& Norris 1997), as well as unpublished radio interferometric 
data showing the nucleus to be unresolved on the 200 pc scale 
(R.R Norris, 2004, private communication). 

The non-detection of several common mid-infrared AGN 
tracers in our Spitzer-l^S high resolution spectrum may hence 
indicate that the AGN coronal line region (CLR) is small 
or absent or, alternatively, that our line of sight to the AGN 
CLR is strongly obscured. The latter would be in agreement 
with the detection of a deep 9.7 ^va silicate absorption fea- 
ture (Spoon et al. 2004, Fig 2), which is unusual for classical 
seyfert galaxies and quasars (Spoon et al. 2007; Hao et al. 
2007). Even stronger obscuration, along all lines of sight, is 
suggested by Lutz et al. (2004), who concluded from a study 
of the 4.7 CO gas absorption band that the strong detection 
of this feature in IRAS F00183-711 1 (Spoon et al. 2004) and 
its non-detection in local AGNs observed with ISO indicates 
that the central power source in IRAS FOO 18 3-7 11 1 must be 
fully covered rather than blocked by an edge-on torus. 

5.2. Excitation conditions in the outflow of 
IRASFOOl 83-7111 

To gain insight into the processes responsible for the pres- 
ence of high velocity gas in IRAS F00183-71 1 1 we compare 
IRAS FOO 18 3-71 11 to the only other two ULIRGs in our 
sample that show pronounced blue wings in their 12.81 /im 
[Ne II] and 15.56 /xm [Ne III] line profiles: IRAS 12127- 
1412NE (Fig. 6) and IRAS 13451H-1232 (Fig. 7). 

Clear evidence for unusual excitation conditions in 



IRAS FOO 18 3-71 1 1 emerge if we compare the R\H emis- 
sion and low-ionization neon line emission for AGNs, 
ULIRGs and starburst galaxies from the various Sp/fzer-IRS 
GTO programs, (U)LIRGs from Spitzer program 1096 and 
IRAS 12127-1412NE. The PAH luminosity should track the 
low-ionization neon line luminosity to within a factor two, if 
predominantly excited in star forming regions. This was in- 
deed found to be the case for the sample of ULIRGs published 
by Farrah et al. (2007). In Fig. 8 we plot the ratio of these two 
quantities (the combined [Ne ll]-i-[Ne III] luminosity and the 
combined luminosity of the 6.2 and 1 1 .2 pm PAH features) as 
a function of increasing apparent depth of the 9.7 pm silicate 
feature (i.e. the negative of the silicate strength, S.,,;; Spoon 
et al. 2007). The symbol size scales linearly with the [Ne 
v]/[Ne II] ratio above a threshold value of 0.4. This distin- 
guishes ULIRGs with a strong AGN contribution to the [Ne 
ll]H-[Ne III] luminosity (e.g. Mrk463; FaiTah et al. 2007) from 
those without. As expected, for the sources lacking a strong 
AGN signature the average value of neon/PAH has a reason- 
ably low scatter of a factor 1.4'^, in agreement with the scal- 
ing relation between neon and PAH flux for ULIRGs found by 
Fan-ah et al. (2007). The notable exception is IRAS F00183- 
7111, which appears to be at least 8 times overluminous in 
the neon lines compared to its combined PAH luminosity. 
Another exception in the Sj,7<-1.5 regime is IRAS 12127- 
1412NE, which is overluminous in neon/PAH by a factor of 
3. At Sj,7>-1.5, the situation is less straightforward. In this 
regime there are several sources with large [Ne v]/[Ne ll] ra- 
tios for which the AGN contributes significantly to the [Ne 
ll]-i-[Ne III] luminosity, thereby raising the neon/PAH ratio 
(Fig. 8). The only ULIRG in this regime which has an elevated 
neon/PAH ratio despite a low [Ne v]/[Ne ll] ratio of 0.28 is 
IRAS 1 345 iH- 1232, at >6 times the average neon/PAH ratio 
(Fig. 8). This brings the number of ULIRGs with unusual ex- 

In this computation we excluded the outliers IRAS FOO 183-71 1 1, 
IRAS 12127-1412NE and IRAS 13451+1232. Keeping them included would 
increase the scatter to a factor 1.7 
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Fig. 9. — Mid-infrared diagnostic diagram for the line ratio 15.56 ^tm [Ne 
III]/ 12.81 /im [Ne II] versus 14.32 /im [Ne V]/ 12.81 ^im [Ne II]. Plotted 
are predicted line ratios as computed by MAPPINGS version Illq (Allen 
et al. 2008) for combined precursor+shock models at solar metallicity at a 
gas density n=100cm"' (upper panel) and n=10cm"^ (lower panel). Dif- 
ferent colors (ranging from purple to red) denote different magnetic field 
strengths, ranging from 10"'^ to 100 fiG. Shock velocities range from 200 to 
1000 km s"' and are labeled along the red curve. Overlaid are the positions 
of IRASF00183-7111, IRAS 12127- 1412NE and IRAS 13451+1232 for the 
integrated fluxes (gray disks) and for just the blue shifted component (black 
disks). Gray squares mark the positions of the nuclear spectra of the classical 
AGNs NGC 1068 and NGC4151 (Sturm et al. 2002) 



citation conditions to three. Since these same three ULIRGs 
are also the only ULIRGs in our sample that show pronounced 
blue wings in their 12.81 /im [Ne ll] and 15.56 /im [Ne III] 
line profiles, it seems plausible that the high neon/PAH ratios 
for these sources could be the result of additional neon line 
emission contributed by these blue wings rather than unusu- 
ally weak PAH emission in these three sources. 

In Fig. 9 we show a mid-infrared diagnostic diagram based 
on the 12.81 ^m [Ne ll], 15.56 ^m [Ne III] and 14.32 ^m 
[Ne v] lines. In the diagram we plot shock excitation mod- 
els from the MAPPINGS Illq library (Allen et al. 2008) 
for solar metallicity gas with densities n=100cm"^ (upper 
panel) and n=10cm~^ (lower panel). The model tracks plot- 
ted include the ionizing effects from UV and soft X-ray 
photons in the shock cooling zone on the gas ahead of the 
shock front. Overplotted are the positions of IRAS F00183- 
71 11, IRAS 12127H-1412NEand IRAS 1 345 1h- 1232 from our 
Spitzer observations. Black filled circles indicate the re- 
sults for just the broad blue-shifted emission line com- 
ponents, the gray versions denote the results based on 
the velocity-integrated line fluxes (i.e. the narrow and 



broad components combined). For IRAS FOO 18 3-71 1 1 and 
IRAS 12127H-1412NE the low observed [Ne lll]/[Ne ll] line 
ratios are consistent with high-velocity shocks (v>500 km 
s"') in high-density environments (n>100cm"-^) with low 
magnetic field strengths (<l/iG). At lower gas densities 
(n=10cm"^), only the highest velocity shocks (v=1000 km 
s"') at the lowest magnetic field strengths (<0.1 /iG) are con- 
sistent with the observed [Ne lll]/[Ne ll] line ratios. Below 
this density the observed line ratios cannot have an origin in 
shock emission. Given the absence of density measurements 
we cannot rule out either scenario. However, as we will ar- 
gue later, if the neon line emission originates at the base of 
the optical outflow, gas densities of n>100cm~^ may not be 
unreasonable. In IRAS 13451-1-1232 it is the high observed 
[Ne v]/[Ne ll] ratio which appears to rule out shocks of any 
velocity and density as the source of excitation for the neon 
lines. For the latter source, exposure of the blue shifted gas 
to direct radiation from the AGN seems a likely explanation 
for the presence of [Ne v] emission at the observed ratio to 
[Ne II]. This explanation is supported by the proximity of the 
blue shifted component of IRAS 1 345 1 -i- 1 232 to the positions 
of the classical AGNs NGC 1068 and NGC 4151 in Fig. 9. 

5.3. Strong obscuration on the outflow in IRAS FOOl 83-71 11 

In Sect. 3.3 we analyzed the velocity structure of 
the 12.81 /im [Ne ll] and 15.56 ^m [Ne III] lines in 
IRAS FOO 183-71 1 1 and compared the results to those of stud- 
ies of its optical and near-infrared forbidden lines. This re- 
vealed that the mid-infrared neon emission lines have strongly 
blue shifted components not detected in both the optical 
5007 A [O III] and neai--infrai-ed 1.644 /xm [Fe ll] fine pro- 
files. Given the mid-infrared spectral evidence for strong ob- 
scuration in the nucleus of IRAS F00183-71 1 1 (Spoon et al. 
2004), the simplest explanation for the absence of strongly 
blue shifted gas components in the optical and near-infrared 
emission lines is extinction. We estimate the required ex- 
tinction by assuming line emission at any wavelength A to 
be effectively suppressed by a specific extinction A(A)>2. 
The level of screen extinction required to obscure the high- 
velocity components of the 5007 A [O III] line then amounts 
to A5007 = Ai/>2. If the same screen is responsible for ob- 
scuring the high-velocity components of the 1 .644 fim [Fe ll] 
line, the required extinction increases to Av>10 (Ai.64/jm>2). 
Then to still be able to see the strongly blue shifted 12.81 /im 
[Ne II] emission through this screen, the maximum extinction 
is constrained to Av<50 (An.%\^j.m<^)- This implies a col- 
umn density of the obscuring screen of Nh = 2-10 x 10^^ 
cm~^, assuming Nh = 1.9 x 10^' x Ay. 

In the absence of spatial information for the mid-infrared 
neon line emission, the location of the screen blocking the 
optical and near-infrared line emission cannot be determined 
with any certainty. However, it is plausible for the screen to 
be situated at the base of the outflow traced in 5007 A [O III]. 
The observed strong decrease in blueshift of the 5007 A [O 
III] line from 3" to 0" East of the nucleus (Fig. 12 of Heck- 
man, Armus & Miley 1990) may hence indicate that the ma- 
terial that obscures the most highly disturbed optical (but not 
mid-infrared) Une emission is located within the inner lOkpc. 

5.4. The evolutionary state of IRAS FOO 183-71 11 

IRAS FOO 18 3-71 1 1 is a clear outlier in a mid-infrared di- 
agnostic diagram which combines the equivalent width of the 
6.2 /.tm PAH emission feature and the silicate strength in the 
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so-called 'fork' diagram (Spoon et al. 2007). In this diagram 
most ULIRGs are found along one of two branches. One 
branch connects starburst-dominated mid-infrared spectra 
(e.g. NGC 7714 and M 82) to spectra of deeply buried central 
sources, characterized by very deep silicate features and the 
(virtual) absence of PAH emission (e.g. IRAS 08572H-3915). 
The other branch connects starburst-dominated sources with 
featureless hot dust-dominated sources, which are mostly 
AGN dominated. 

IRAS FOG 183-71 11 is located in the scarcely populated 
region in between the tips of the two branches and has 
IRAS 12127H-1412NE as its nearest neighbor. Their location, 
away from other deeply obscured sources, either implies a 
more than 10 times lower 6.2 fim PAH equivalent width than 
other ULIRGs with similar silicate strength, or a less pro- 
nounced silicate feature than other deeply obscured ULIRGs. 
Spoon et al. (2007) have proposed that sources in between 
the two branches may be AGN evolving from a fully covered 
obscuring geometry to an incomplete or clumpy obscuring ge- 
ometry, with a decrease in apparent depth of the silicate fea- 
ture as a result. In this scenario, the detection of strongly blue 
shifted neon line emission would then trace a unique phase 
in the (partial) disruption of the obscuring medium around 
the powerful AGN. Alternatively, the isolated positions of 
IRAS F00183-7111 andlRAS 12127H-1412NEin thefork di- 
agram may be the result of unusual lines of sight into these 
systems. The unusual mid-infrared excitation conditions and 
kinematics may then be associated to that. 

If sources with mid-infrared kinematic properties and 
neon/PAH ratios similar to IRAS FOO 183-71 11 and 
IRAS 12127-H1412NE are indeed AGN evolving from 
full to incomplete or clumpy obscuration, there should be 
observational characteristics linking them. However, the 
only ULIRG in our sample with these same mid-infrared 
properties, IRAS 13451H-1232, has a mid-infrared continuum 
spectrum which is very different from the other two sources, 
displaying only a weak silicate absorption feature on an oth- 
erwise featureless power law spectrum. This places it firmly 
on the AGN branch in the fork diagram, far off from the other 
two sources. Evidence supporting an unusual evolutionary 
stage is presented by Holt, Tadhunter & Morganti (2003). 
Based on the shape of its radio SED and a kinematic study 
of its optical emission line spectrum, they conclude that 
IRAS 13451-H1232 must be a young compact radio source 
whose radio jets are expanding through a dense nuclear dust 
shell, entraining gas and dust from this shell. They infer the 
reddening to be highest toward the broadest (most disturbed) 
most blue shifted kinematic component in their optical line 
study, reminiscent of our findings for IRAS FOO 183-71 11 
using optical, near and mid-infrared line profiles. 

6. CONCLUSIONS 

We have discovered high-velocity ionized gas with veloc- 
ities ranging from -3500 to -h3000 km s~^ as part of a nu- 



clear outflow from the ULIRG IRAS F00183-7111. The gas 
is traced by the mid-infrared 12.81 fim [Ne II| and 15.51 fim 
[Ne III] lines, which are dominated by blue shifted compo- 
nents. Optical and near-infrared spectroscopic studies show 
no evidence for similar high- velocity gas components in for- 
bidden lines at shorter wavelengths. We interpret this as the 
result of strong extinction (Av=10-50; Ai2.8i/im<l) on the 
blue shifted gas, identifying the base of the optical outflow 
(<3" East) as its most likely origin. 

Unusual excitation conditions in IRAS FOO 183-71 11 are 
further implied by a comparison of the combined [Ne ii] 
and [Ne III] luminosity and the combined 6.2 and 11.2/im 
PAH luminosity for a large sample of AGNs, ULIRGs and 
starburst galaxies. This reveals the neon line emission in 
IRAS FOO 183-71 11 to be at least 8 times stronger compared 
to its PAH emission than the typical ratio for the sample. 

Two other ULIRGs in our sample of 56 display 
mid-infrared excitation and kinematic properties sim- 
ilar to IRAS FOO 183-71 11: IRAS 12127-1412NE and 
IRAS 1345 lH-1232. Both sources have elevated neon/PAH ra- 
tios and exhibit pronounced blue wings in their 15.56 fim |Ne 
III] line profiles. Wings are also present in their 12.81 fim [Ne 
II] profiles, but are less pronounced, indicating blue shifted 
gas to have a higher [Ne lll]/[Ne ll] ratio than gas closer to 
systemic velocity. IRAS 13451-1-1232 even shows strongly 
blue shifted 14.32 fim [Ne v] emission. 

While for IRAS 13451-1-1232 the observed [Ne lll]/[Ne ll] 
and [Ne v]/[Ne II] ratios indicate exposure of the blue shifted 
gas to direct radiation from the AGN, for IRAS F00183-7 1 1 1 
and IRAS 12127-1412NE the observed ratios are consistent 
with an origin in fast shocks (v>500 km s~^) in an environ- 
ment with gas densities n>10cm"^. 

IRASF00183-7111 and IRAS 12127-1412NE are clear 
outliers in the mid-infrared diagnostic diagram of silicate 
strength and 6.2 /xm PAH equivalent width, the so-called 
'fork' diagram (Spoon et al. 2007). The sources are found 
in between the two branches which separate deeply buried 
nuclei and sources with incomplete or clumpy obscuring ge- 
ometries. Spoon et al. (2007) have proposed that ULIRGs in 
between the tips of these branches may be AGN in the pro- 
cess of disrupting their fully covered geometry before settling 
on the AGN branch. IRAS 1 345 1h- 1232 would then already 
be there, given its nearly featureless mid-infrared power law 
spectrum. 

The results from our mid-infrared kinematic and excitation 
study are consistent with the above evolutionary scenario in 
that the strongly blue shifted [Ne ll] and [Ne III] emission 
(from shocks) may trace the (partial) disruption of the obscur- 
ing medium around buried AGNs in IRAS F00183-71 1 1 and 
IRAS 12127-1412NE. The detection of strongly blue shifted 
[Ne V] emission in IRAS 13451-H1232 may then indicate the 
lifting of the obscuration around the AGN coronal line region 
to be furflier advanced in this ULIRG than in IRAS FOO 183- 
71 1 1 and IRAS 12127-1412NE, where this line is undetected. 



APPENDIX 
Molecular hydrogen lines 

Our measurements of the S(l), S(2), S(3) and S(4) pure rotational molecular hydrogen lines allow us to investigate the prop- 
erties of warm molecular hydrogen in IRAS FOO 183-71 1 1. Unfortunately, a measurement of the important S(0) line is lacking 
as it is outside the IRS-LH bandpass. In Fig. 1 we show the excitation diagram for molecular hydrogen for the four measured 
lines. Given the placement of the S(3) line at the peak of the silicate absorption feature, we can attempt an extinction correction 
on the four data points under the assumption that a) the excitation temperature should be constant or increase monotonically with 
increasing upper level energy Ej,(i/,J), b) that the ortho-to-para ratio is 3 (i.e. assuming LTE), and c) that the cooler gas probed by 
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Fig. 1. — Excitation diagram for warm molecular hydrogen in IRASF00183-7111 based on the line fluxes for the pure-rotational transitions S(l)-S(4). 
Different shades of gray are used to distinguish different levels of dereddening: black: A\-25, dark gray: Av=10, light gray: Av=0. Error bars are too small to 
be visible 



the S(l) and S(2) lines is more or less co-spatial with the warmer gas probed by the S(3) and S(4) lines. Adopting the local ISM 
extinction law of Chiar & Tielens (2006), we then find that the extinction on the H2 lines ranges between Av=15 and 25. This is 
far lower than the extinction on the continuum source, which is Av>5A based on the apparent depth of the 9.7 jim silicate feature 
(t9.7=3.0; Spoon et al. 2006) and Av/r9.7=18 (Roche & Aitken 1984). Based on the two detected ortho-Hi lines, it is possible to 
infer an excitation temperature for the warm molecular hydrogen. The value ranges from 330 to 390 K for extinction corrections 
ranging from Ay=0 to 25. Assuming all H2 levels to be in LTE at 330-390 K, we derive a total warm molecular hydrogen mass 
of 8-20 X 10** M0. This value becomes a lower limit if the excitation temperature based on the upper level populations of the 
S(0) and S(l) lines is lower than those based on the S(l) and S(3) lines - which is usually the case (e.g. Higdon et al. 2006). 
Unfortunately, we cannot determine the warm to total molecular hydrogen fraction and compare this to other ULIRGs, as no 
millimeter CO line measurement exists for IRAS F00183-71 1 1 from which to infer the cold molecular hydrogen mass. 

Non-detection of absorption bands of C2H2 and HCN gas 

Several ULIRGs show absorption bands of C2H2 and HCN gas in their 13.7-14.1 /^m Spitzer IRS-SH spectra (Spoon et 
al. 2005; Lahuis et al. 2007). Detections are especially strong among ULIRGs with pronounced silicate absorption features, 
like IRAS 20100-4156 and IRAS 15250H-3609. Far weaker detections are reported for ULIRGs with AGN-like spectra, like 
IRAS 05189-2524, IRAS 01003-2238 and Mrk231 (Lahuis et al. 2007). Detections at these levels require spectra to have signal- 
to-noise ratios of at least 100 in their continua. The IRS-SH spectrum of IRAS F00183-711 1 does not show any indication for 
the presence of either feature, despite a continuum S/N of 70. This ttanslates into upper limits on the C2H2 and HCN gas colunan 
densities of N(C2H2)<1.2x lO^^cm'^ and N(HCN)<3.0x lO^^cm'^ for an adopted excitation temperature of 300K (F. Lahuis, 
private conomurucation). 
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